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Within a quasipotential framework a relativistic analysis is presented of the deuteron current. 
Assuming that the singularities from the nucleon propagators are important, a so-called equal time 
approximation of the current is constructed. This is applied to both elastic and inelastic electron 
scattering. As dynamical model the relativistic one boson exchange model is used. Reasonable 
agreement is found with a previous relativistic calculation of the elastic electromagnetic form factors 
of the deuteron. For the unpolarized inelastic electron scattering effects of final state interactions 
and relativistic corrections to the structure functions are considered in the impulse approximation. 
Two specific kinematic situations are studied as examples. 
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I. INTRODUCTION 



Electron scattering processes have been crucial in extracting detailed information about the nuclear interaction. 
Various exclusive Jl],0,[|] and inclusive |^,^| breakup reactions for the few nucleon systems have recently for example 
been performed at moderate momentum transfer, where separations of the longitudinal and transverse response 
functions have successfully been achieved. The deuteron, as the most simple nuclear system, is of special interest 
because exact calculations are in principle feasible for this case. Such a system may serve as testing ground for 
theoretical models describing the nuclear dynamics and the e.m. operators used to study these reaction processes. In 
most of these studies it is implicitly assumed that the constituents of the system behave nonrelativistically and that 
the e.m. interaction can be treated in an independent way. With increasing momentum and energy transfer effects of 
relativity are expected to play an important role and a relativistic description should be essential in these kinematic 
regions. Moreover, in view of gauge invariance a consistent treatment of the e.m. interaction is called for. 

In this paper we discuss a framework which is well suited for a relativistic analysis of both the clastic and inelastic 
electron scattering off the deuteron. It is based on the relativistically covariant field theoretical Bethe-Salpeter equation 
approach. Gauge invariance is satisfied through the Ward-Takahashi identity. Assuming a theoretical description in 
terms of nucleonic and mesonic degrees of freedom the elastic e.m. properties have already been studied. Within an 
one boson exchange (OBE) model B a consistent relativistic treatment of both the e.m. current and the nucleon- 
nucleon interaction can be realized 0, including also meson exchange current( MEC) contributions from the pn^f and 
loe^i currents In the actual calculations a relativistic quasipotential approximation of the Blankenbecler-Sugar- 
Logunov-Tavkhclidzc (BSLT) has been used. 

The disintegration experiments of the deuteron yields additional information as compared to the elastic e.m. scat- 
tering process. A larger kinematic region can be tested and the theoretical models are more complicated due to final 
state interaction (FSI) and other MEC contributions as is the case in elastic scattering. Another point of interest is 
that, depending on the type of experiment, the cross section can be separated in more than 2 structure functions. 
The inclusive experiment is described by the longitudinal Rl and transverse Rt structure functions. For the exclu- 
sive experiment studied here an additional structure function Rtl is measured, which is an interference between the 
longitudinal and transverse components in the current operator. Because this structure function depends on other 
reaction processes as Rl and Rt, new information can be obtained. 

In our calculations complete knowledge of the FSI and the deuteron wave function is required. As mentioned 
earlier we use the OBE model of Ref. || . Although the complete model can be studied we use here a quasipotential 
approximation to it, where both particles are treated in a symmetrical way. This BSLT approximation p|,|lO|| was also 
used in the model of the elastic form factors. Here we'll use the same approximation for the initial and final state, 
but use a different approximation for the nucleon propagators occurring in the e.m. vertex. We assume that only 
a relative energy dependence is in these propagators and not in the two-nucleon vertex functions. As a result the 
relative energy dependence can explicitly be integrated out. Because the BSLT choice is used for the nucleon states 
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this is an equal time (ET) approximation to only the current operator. To examine the sensitivity of the predictions 
on such a ET prescription for the e.m. current matrix elements we calculate the e.m. form factors in the elastic case. 
Comparing to earlier calculations || we find only small differences. The ET choice has the advantage that it can 
also be extended in a systematic way to describe the case of inelastic electron scattering. In the full Bethe-Salpeter 
theory a conserved deuteron current can be constructed (see Refs. 0,|ll|), at least for the case of strong form factors 
at the meson-nucleon vertices, which depend only on the momentum of the meson. The same can be achieved in 
a systematic way in the ET approximation, i.e. MEC contributions can be constructed for the breakup reaction to 
ensure current conservation. 

The paper is organized as follows. In the next two sections we summarize the relevant expressions, characterizing 
the inelastic electron scattering processes and describe the dynamical model used in this work. In particular, we 
present a new fit in this model to the experimental phase shifts, to discuss the effect of including the negative energy 
states. In sections IV and V the relativistic impulse approximation and the elastic e.m. form factor calculations 
are described in the ET approximation. Section VI deals with the formalism of the e.m. breakup of the deuteron 
in this approximation. In all these calculations a fully relativistic e.m. current operator with on shell form factors 
is used. In section VII two kinematic situations are studied using this relativistic formalism. For the conventional 
nonrelativistic (NR) models different forms of e.m. operators have been constructed in the literature. These models 
differ predominantly in the use of F\ or Ge for the e.m. nucleon form factors and the differences are obviously of 
relativistic order. We compare these predictions with our fully relativistic analysis and show that in the kinematics 
studied the predictions can be very different. Most of the differences can be traced back to the choice of the e.m. 
operator. Although the experiments are at moderate momentum transfer, large relativistic effects are in particular 
found in the interference structure function Rtl, showing the failure of the NR models discussed here. The next 
leading order correction to these operators is given and shown to represent reasonably well the relativistic predictions, 
at least in the kinematics considered. In the kinematic situations discussed in this section the MEC contributions 
which are needed to ensure gauge invariance are of minor importance. In a next paper we will discuss in detail our 
MEC calculations in this relativistic model. 



II. DEUTERON STRUCTURE FUNCTIONS 



In this section we summarize some useful formulae describing the electrodisintegration process of the deuteron. The 



conventions of Bjorken and Drell |12| are used. We confine ourselves in this paper to the description of unpolarized 
electron scattering on an unpolarized target and follow closely the work of Donnelly and Raskin jL3|. The analysis 
can readily be extended to the case of polarized scattering. 

For definiteness, let us consider the breakup of a deuteron with total momentum P by an electron into a free 
neutron-proton (np) pair, characterized by the four momenta p n — (p n ,E n ) and p p = (p p ,E p ) respectively. In 
the one-photon exchange approximation, the differential cross section in the lab system, being the rest frame of the 
deuteron i.e. P = (0, Mo), can be written as 

d ° - c ^ E ^-0f d3 p? d3 p^ k +p-x- p') e ^\f^\ 2 , (2.i) 

if 

with C np — M%f/(E n E p ) and where is the electron and J M =< f \ j£ \ i > the deuteron electromagnetic current 
matrix element between the initial and final state of the nuclear system. The total four momentum of the np pair 
is assumed to be given by P', while q is the four momentum of the virtual photon satisfying q^ — fc M — k'^ = (q, u>), 
hp = (k, e) and k'^ = (k', e') being the four momenta of the incoming and scattered electron. Momentum conservation 
at the photon-deuteron vertex gives P + q = P' . The normalization of the outgoing np pair state is such that when 
the final state interaction is neglected it is given by < p pi p n \ f >~ v>(p p )u(p n ), u being Bjorken-Drell spinors. In 

Eq. ( pi] ) ^2 indicates averaging and summing over both the electron and nuclear polarizations in the initial and final 
if _ 

state respectively. An electron and hadronic current tensor 

»w = E tifi> w ^ = c ™p E J t J - ( 2 - 2 ) 

if if 

can now be introduced such that |j^ M | 2 = rj^W^. The expression for 7/ M „ is well known jl2| and is given in the 

if 
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ultra relativistic limit ( m e w 0) by 



V = ^(k^K + k^-g^k-k'). (2.3) 



Introducing the relative momentum p' = (p p — p n )/2 of the np pair and integrating over P' yields for the cross section 



ml e 4 d 3 k'd 3 p', 



At this point the concept of the structure functions can be introduced. We may choose a coordinate system, where q 
is along the z-axis, while the y-axis is defined along k x k'. Let us now define in this lab frame J± = -F^j (Jx i iJy)i 
and use current conservation q ■ J = to e limi nate the longitudinal third component J3 of the current in favor of the 
charge J3 = u/qJo. With the help of Eq. (2.3) the contraction can be written in the form 

"H^W^ = v a (v l R l + v t Rt + v tt Rtt + vtlRtl) , (2.5) 

with vo = 4ee'cos 2 \Q e , 8 e being the scattered electron angle. Moreover, the electron-kinematic factors are given by 

(2 \ 2 2 2 2 1 2 

^ J , v T — ~\ < ^i + tan 2 |0 e , vtt = and vtl — "^f^y + tan 2 \0e- The nuclear structure functions 

R a can be expressed in the matrix elements of the various current components as fl3| 

Rl = C np \ Jo| 2 , 
R T = C np (\J + \ 2 + \J-\ 2 ), 
Rtt = 2C np Re (j; J J) , 

R TL = -2C np Re [J * ( J+ - J_)] . (2.6) 

It should be noted that C np can be absorbed in the Bjorken-Drell spinors occurring in the e.m. current matrix 
elements. The resulting spinors become essentially those defined by Kubis [Q, described in Appendix A. In the 
remainder of the paper these spinors will be used. Introducing the Mott cross section OMott = ^t—vq, with the fine 
structure constant given by a = e /(4tt), the cross section can be written as 

da = a mm {vlRl + v T R T + v TT Rtt + vtlRtl) S(M d + lu - E p - E n )de'dn' e d 3 p . (2.7) 

This is still a general expres sion , except that current conservation of the hadronic current is assumed. 

All the quantities in Eq. (|]?]) are determined in the lab frame. Considering the total four momentum of the np 
pair P' = (q,E np ) (with q = p„ + p p and E np = E n + E p ) an invariant mass M np of the pair can be defined as 

E np — ^/q 2 + M 2 p . Since the final state interaction of the final np pair can most readily be determined in its cm. 

system (p' cm = (0,M„ p )), we prefer to express Eq. (2.7) in terms of the structure functions R^ 11 , evaluated in this 
frame. To make a distinction between lab and cm. variables we use the convention that all variables in the cm. 
frame are labeled by cm. The relation can easily be found from 

j; m = £ l rj v . (2.8) 

where C" describes the boost between the lab and cm. frame. An explicit expression for the Lorentz transformation 
can be found using the relation P p m = C^P' V with 

VT+V -y/rj 
10 
^ I 1 



Ai I n n 1 n I ; (^-9) 



where we have 



-07 VT+rj 



Vv = ^,ViT-v = ^. (2.10) 

lvl np lvl np 



Since q is along the z-axis, the transverse components of the deuteron current are not affected by the Lorentz boost. 
For the charge component we get J§ m = y/1 + t/Jq — ^/rjJ^. Together with current conservation q - J = q cm ■ J cm = 0, 
we find a direct relation between the charge components of the current in the two frames 
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Jo = A 4r J ° cro ^ 

Mo 

We now turn to derive the differential cross section describing the exclusive die, e'p)n and inclusive (e, e') reactions. 
In view of the Bjorken and Drell covariant normalization of the free np spinors, the transfor mati on property of the 
C np coefficient and the invariant volume element d 3 p' cm / E^V 1 = d 3 p' /E p > we may rewrite Eq. (2.7) as 

da = ouott (vlRl + v t Rt + v TT R TT + v TL R TL ) 8(EfF + co cm - 2Ep n )de'd{l' e d 3 p' cm , (2.12) 
where the structure functions can be expressed in the corresponding observables in the cm. system 

r> f Mnp \ r?cm p rtcm r> _ rycm p ( Mnp \ ncm fry in\ 

The exclusive reaction d(e, e'p)n is represented by a five fold differential cross section. In view of the invariant volume 
element we immediately find from Eq. (2.12) 

d 5 a ( \ 

CMott [vlRl + v T R T + vttRtt + vtlRtl I , (2-14) 



de'dSl'dQ' 



where we have absorbed the Jacobian J into the definitions of the exclusive structure functions R a = JR a . For the 
Jacobian we have J = d " n T \p' cm E'5 n with 

dQ ' Cm - P> dE * - P ' 1 d-^St***)- 1 (2.15) 



dfl' p' cm dE c v V l p' cm Vl+T] V V 1 + J 1 P' 

Hence, in principle the unpolarized exclusive reaction can be calculated by determining the four structure functions. 
For the (e, e') process the integration over the momenta of the final np pair has to be performed. This is actually 
done in the cm. frame of the np pair. Integration over the relative momentum p' cm yields the inclusive (e, e') cross 
section (in the lab frame) 

d 3 <r , — — N 

^7^y = VMott [VlRl + V T Rt) , ( 2 -16) 

where the inclusive structure functions R a arc defined by 

f 1 1 

R a = 2nJ R a -p' cm E^ n d cos 0' cm . (2.17) 



The momentum p'™ is determined by the energy conserving S function in Eq. (2.12), i.e. from Ey m = ^ (E^ 1 + uj cm ). 
The contributions to the inclusive cross section from Rtt and Rtl, which are proportional to cos2cf>' cm and cos(j)' cm 
respectively, vanish due to the 0' cm -integration. In section VI we will discuss this <fi' cm dependence in more detail. 



III. THE FORCE MODEL 



In determining the nuclear structure functions we in principle need to know the deuteron vertex function and the 
half off-shell nucleon-nucleon (NN) t-matrix. The nuclear interaction used by us is based on the one-boson-exchange 
(OBE) model. Following Fleischer and Tjon Q it is assumed to be described by the exchange of 7r, p, cj, 77, e (or a) and 
5 mesons. A strong meson-nucleon form factor of the monopole type F{p 2 ) — ^t_ p i is used to regularize the large 
momentum behavior. Within the relativistic field theory two-particle scattering can be described by the scattering 
t-matrix <p, which satisfies the Bcthc-Salpctcr equation. The inhomogencous Bethe-Salpeter equation for the t-matrix 
has the form 



<t>(p',p;P) = V(p',p)-^- i J d 4 k<j ) (p\k;P)S 2 (k,P)V(k,p) 



V(p',p)-^ I d 4 kV(p',k)S 2 (k 1 P)^k lP ;P), (3.1) 
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where 



{\ff + 1 f,+ M N )V>{\l? -i> + M N )&>Go, (3.2) 



with G = [(iP + p) 2 -A/2 r + ie]- 1 [(iP-p) 2 -M2 r + ie]~ 1 . In Eq. Q S^(p,P) is the free propagator of the n-th 

nucleon. The t-matrix also determines the deuteron vertex function $^f f ' > (p,P) (where M is the polarization of the 
deuteron), corresponding to the residue at the deuteron bound state pole at P 2 = Mp. For P 2 « Mq we have 

0(p', P! P)=£ 02 _m2 + regular terms. (3.3) 

M D 

Alternatively, $d satisfies the homogeneous equation 

^d\p,P) = ~ J d i kV(p,k)S 2 (k,P)^ I \k,P). (3.4) 

Although a field theoretical Bethe-Salpeter analysis is in principle possible for the case that V is given by the OBE 
model, calculations are highly non trivial because of its analytic structure and computational complexity. To simplify 
this we use a relativistic quasipotential approximation. In the quasipotential framework the two-particle propagator 
is replaced by one where the relative energy variable is constrained (for a review see (l5[ ). Here we use the choice 
which treats the two nucleons in a symmetric way pUlpfl . It is given by 

Go -> inS( Pa )Gi SLT = nrS( Po ) — \ - — * -y 2 , (3.5) 

where E = \Pq and E p — M'fj + p 2 . In this approximation the full two-particle propagator S2 SLT , containing the 
spinor structure (see Appendix A), is given by 

Si SLT = {\& + i + M N )V{\P - i + M N )^Gi SLT = (E p - E)SM (p, P)SW (p, P). (3.6) 



In this approximation the inhomogeneous Bethe-Salpeter equation ( p.l| ) simplifies to 

0(p,p';P) = V(pJ) + d 3 k<j ) (p,k;P)S 2 BSLT (k,P)V(k,p'), (3.7) 

where p,p' and are the four vectors p,p' and fc, subject to the condition that their fourth component vanishes 
in the two-nucleon cm. frame. To reconstruct the off-shell scattering wave function needed in the study of FSI 
effects a partial wave representation is used for the full NN amplitude. Following Ref. ||, the helicity basis is used. 
The representation is briefly di scus sed in Appendix A. The on-shell amplitude is simply obtained by taking p = p' 
with a/p 2 " = \/s = 2E P . Eq. (3/7) can be sol ved in a partial wave analysis. For details we refer to [||. Because 



the quasipotential approximation is used, Eq. (3.7) reduces essentially to a coupled set of one-dimensional integral 
equations. Besides the physical (+, +) positive energy states, also (— , — ) states and even and odd combinations of 
the (+, — ) energy states occur Q . 

Starting from the original fit of Ref. || (set A) for only positive energy spinor states, the meson-nucleon coupling 
constants g e and gs were adapted to the case that all the spinor states are included to get a reasonable fit to the 
experimental phase shifts of Arndt et al |Tq] . The resulting phase shifts up to J = 2 are shown in Fig. [| (set B), where 
the spectroscopic notation 2S+1 Lj has been used. For comparison the results of fit A are also plotted. To see what 
the effects are of the negative energy spinor states, we have switched these states off, keeping the same strength for 
the coupling constants. The results are shown as the dotted-dashed lines. Finally, an attempt was made to vary the 
meson coupling constants for the case of only positive energy states to reproduce the phase shifts as obtained with all 
spinor states included. We indeed obtain phase shifts (set C) which are not distinguishable from the set B in fig. |]. 
The sets of coupling parameters are given in Table I. 
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IV. IMPULSE APPROXIMATION 



In order to construct a quasipotential approximation for inelastic electron scattering we start with the current 
operator as given in the Bethe-Salpeter formalism. The contributions to the current operator in the impulse approx- 
imation (IA) are given by the Feynman graphs shown in Fig. |^ and Fig. ^ in the case of inelastic electron scattering. 
In the contributions (a) and (b) ((c) and (d)) in Fig. || the photon couples directly to the proton (neutron) which is 
knocked out without any interaction with the other nucleon. These contribution are called the plane wave impulse 
approximation (PWIA) and the Born term respectively. The set of graphs in Fig. [| corresponds to the contributions, 
where the outgoing nucleons interact after the virtual photon has been absorbed. These are the so-called final state 
interaction (FSI) contributions. In our calculations we will absorp the Born contribution in the FSI contributions. As 
is seen from the figure the FSI is expressed in terms of the half-off-shell NN t-matrix cf>. It satisfies the inhomogeneous 



equation given in Eq. (3.1). Including the Fermi character of the nucleons, the free np pair with total momentum P' 



and relative momentum pj is given by the antisymmetric combination of two free Dirac spinors with helicities X n 

< Pf ,P'\=u Xl (±P'+p f )ux 2 (\f"-Pf)A, (4.1) 

where A is the antisymmetrizer. As mentioned before we use the Kubis spinors. Since the outgoing particles are 
on-shell, we have to satisfy the two conditions 

Pr P' = 0, = E iP , +pf + £ 4p ,_ p/ . (4.2) 

The PWIA and FSI contributions can be combined by introducing the scattering wave function for the two nucleons 

* np (p f ,p';P') =< Pf ,P'\ [4^<SV -Pf) + k P (Pf,p';P')S 2 { P ',P% (4.3) 



where 4> np satisfies Eq. ( |3.1| ), subject to the condition that both outgoing particles are on mass shell. As a result, the 
current operator in the IA as represented in Figs. ^ and ^| can be written in the form 

4 A = *1 Md E / d^ np {jp } ,p'-,P')Tf{ q )S^{p^,P)^ D {p^,P), (4.4) 

where p^ = p' — \<l and p^ = p' + ^q. Because of four momentum conservation we have P 1 = P + q. For the jNN 
vertex for the k-th nucleon we take 

r£>(g) = F[ k \ q )^ + J^F« t \ q )*Wq-> > (4.5) 

with F^ k) = F.s + T^ k) FX. In the actual studies we assume that the e.m. nucleon form factors F n can be described by 
their on-shell form. To obtain a similar form for the current matrix element as in the elastic case we may introduce 
the corresponding scattering vertex function 

$ np (p,p';P) = ^ np (p,p';P)S^(p',P). (4.6) 



From studies of the two- nucleon functions <£>„ p and $d 17 1 a rather smooth behavior is found in the relative energy 



Po variable. Therefore we assume that an expansion of the vertex functions in Eq. (4.4) around their relative variable 
Po = point is reasonable in practical calculations. Keeping only the lowest order contribution and assuming that 
the resulting $'s can be taken to be the BSLT vertex function, we get for the e.m. current operator 

4 A = T^ftsr E / dV^ T (p';P05 2 (p^P0^«5W(p«,p)$^(pW;P). (4.7) 

^ ' i—l 

Since the only po dependence is in the nucleon propagators we may carry out the po integration analytically. In 
Appendix B the explicit formulae of the propagator structure are given. The above approximation is essentially an 
equal time (ET) approximation to the current operator, i.e. with zero relative time. 

Let us consider the covariance aspects of the quasipotential description. Under a given Lorentz transformation C 
we have 



G 



V(k,p) = A{£)V{£- 1 k,£~ 1 P )A- 1 {£), 

S w (p, P) = A« (C)S W (£- x P , £~ 1 P)A^~ 1 (£), 
<5> D {p,P)=A{£)$ D {C- l p,£- 1 P), 
(k,p, P) = AOC^pOC- 1 *, £^p, £- 1 P)A-\£), 



(4.8) 



with A(£) = A (1 )(£)A( 2 )(£), where A«(£) is the corresponding boost operator in the Dirac space of the i-th nucleon. 
Because the BSLT deuteron and continuum wave functions transform in a similar way as the Bethe-Salpeter ones, 
the transformation property of the deuteron current for the quasipotential case is the same as that of the relativistic 
field theory case. Since 



A( 1 )" 1 (£)rW(g)A( 1 )(£) = C?T<?\£-\), 



(4.9) 



we may conclude that the deuteron current given by Eq. ( [4.7D indeed transform as a four vector. 

We close this section with some remarks on current conservation. As noted previously, for elastic electron deuteron 
scattering the current operator as defined in the Bethe-Salpeter formalism is conserved, provided the kernel of the 
Bcthc-Salpctcr equation satisfies a local property. One of the crucial ingredients in showing this is that the jNN 
vertex satisfies the Ward-Takahashi identity 



q.T® =F$ i) (S® \p^\P')~ S* w \p,P)). 



(4.10) 



The analysis can readily be e xten ded to the case of inelastic electron scattering. Note that in view of Eq. (3T) the 
scattering wave function Eq. (4.3) satisfies the homogeneous equation 



(4.11) 



Substituting the Ward-Takahashi identity in Eq. (4.4) and using Eqs. ( p.4[ ) and (4.11) leads in the same way as in the 
clastic case to 



■r-r 1 = {27r y MD E J d " k j d 4 p^ np {k,P') [v{k,p*),F^\ S 2 (p,P)<S> D ( P ,P). 



(4.12) 



The main difference with the case of elastic electron scattering is that the final two-nucleon state is different from the 
initial state and that the e.m. nucleon form factors in the commutator also contain in this case in addition to the 
isoscalar also an isovector part. If the final state is an / = state the form factors are isoscalar and the commutator 
vanishes as for elastic scattering. In breakup reactions however also the 7=1 channel is present in the final state, 
yielding a non vanishing isovector contribution. Consequently MEC contributions should be added to Eq. (4.12) to 
make it divergenceless. These will be considered in detail within the quasipotential approach in a forthcoming paper. 



V. ELASTIC SCATTERING IN ET APPROXIMATION 



The elastic deuteron current has precisely the same form as Eq. (|17|) except that $ — > $£> . It is worth noting that 
the deuteron current for elastic scattering is conserved in the ET approximation at the level of posi tive energy states. 
This can be seen as follows. Keeping only the positive energy state contributions we get from Eq. (4.7) 



q-J 



I A 



6&TT 5 M D 
2 

X 

i=l 



d 3 k / d 3 p$Z SLT (k;P')Si SLT (k,P' 



[v(k,p%F^] Si SLT (p,P)^ D SLT (p;P), 



(5.1) 



since l/iir J dpoSzip, P) is equal to Sjf T (4p, P) for positive energy states. Because now only the isoscalar parts of 



(%) 

the form factors F± contribute to the current ma trix elements for elastic scattering and these parts commute with 



the potential V, we may indeed conclude that Eq. (|5.1|) is divergenceless. 



7 



Since the normalization condition of the deuteron wave function is related to current conservation, we should expect 
the correct normalization only for positive energy states. The normalization condition for the BSLT vertex function 
is given by 

L_ f rfuZW/u. m ( JL G BSLT (1 „ u\\ 



P m 6m,m> = (2^)3 7 d 6 k^"(k;P) ^ — S^(k,P) ) Qtf >(k;P). (5.2) 



p 2 =m5, 



This can readily be verified to hold. Using the identity 



d -S 2 (p,P) = hsW(p,PhPs 2 (p,P)+sW(p,P) 1 Ws 2 (p,P)}, (5.3) 



dp,,. ^ ' 2' 



we get for the ET current matrix element Eq. (4.7) in the limit q — > 



4 A = j^fr^ J d^ B D SLT (p, P)£rS 2 (p, m B D SLT (P, P), (5.4) 



which is indeed identical to Eq. ( |5.2| ) provided that only the positive energy states are kept. It should be noted that 
current con servation for the more general case can in principle be restored by adding an effective two-body current 
to Eq. ( |5.l| ) to cancel the four divergence of the IA current matrix element. Because the violation is at the level of 
the negative energy state contributions, the correction is of relativistic order. 



The deuteron current (4.7) as discussed in the previous section is obviously different from the version studied in 
||, where the BSLT two-body propagators were used for both the initial and final state propagators. Here we have 
kept the propagator structure and the po-integration is performed only over the propagators, since we have implicitly 
assumed that the po-dependence may be neglected in the vertex functions. We now compare the elastic electron 
scattering results as obtained with the two approaches. The calculations in the ET approach proceed in exactly the 
same way as in |^|,[l8). Using the helicity framework (see Appendix A) we get for the current matrix element 

< P\M'\4 A \P,M > = {2 ^7 Md J dk ° J dkk " J dV k 

x E E E^^^^a^w^W' 

nri' 77i, m' p,p' 

(5.5) 



where we use the notation m = {Ai,A 2 ,p} rh = {Ai, X 2 , p}, with Ai being the helicity of the i-th nucleon and p 
the quantum numbers of the rho-spin of the two-nucleon state. The angular quantum numbers are characterized by 
n = {J,M,L,S}. 

The current matrix elements are evaluated in the Breit frame. Both initial and final states are boosted to their cm. 



frame. On the p-spin basis |++>, | >, |e > and \o > both propagators in Eq. (5.5) can be written in matrix 

form 



S++ 
5__0 
S ee S 

ec 

S oe S oc 



(5.6) 



For the single nucleon propagator we have S++ = S^(k'), S— = S'i 1 ^'), S ee = S ao = (S { +\k r ) + S i } ) (k'))/2 and 
S eo — S oe — (S+ (k 1 ) — S^\k'))/2. The two-nucleon Green function 52- (k) can also be written in the form (5.6) 



with S++ = 5 2++ (fc), 5__ = S 2 __(k), See = 5 OG = (5 2+ _(fc) +5 2 _ + (fc))/2 and 5 eo = 5 oe = (5 2+ _(fc) - S 2 _ + (k))/2. 
For the explicit form of these propagators we refer to Appendix A. We note that the matrices are no longer diagonal 
as is the case for the choice made in j|] where the BSLT propagators are used. 

In Fig. H are shown the results for the electric form factors Fc and Fq. The dashed line and the dotted-dashed 
line give respectively the results when all components or only positive energy states in the deuteron vertex function 
are included. We see that only at high momentum transfer the negative energy components contribute noticeably 
to the form factors. For comparison also are shown the results for the BSLT reduction as discussed in ||. In both 
calculations the Hohler et al. jNN form factors have been used. We see that the dip in the charge form factor 
Fc is sensitive to the model used. At low momentum transfer the models give results in good agreement with each 
other. Only at high q 2 the results differ. In Fig. || are displayed the results for the magnetic form factor Fm- Negative 
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energy components in the wave function yield large contributions to the form factor in the dip region. When these 
components are included, the position of the dip is shifted to higher momentum transfer. The two quasipotential 
reductions give different results when only positive energy states are taken into account. However, inclusion of the 
negative energy states results in a good agreement between the two approximations. The results for the electric form 
factor A and the magnetic form factor B are shown in Fig. ^ using the ET approximation. In the calculations of the 
(p7T7 and ue^) MEC the BSLT propagator has been used instead of the ET one. They only differ (see Appendix A) 
for negative energy states. Calculations in various kinematic situations show that such a replacement doesnot affect 
the MEC results noticeably. The complete ET calculations of A yield essentially the same results as those with the 
BSLT approximation, except for some small differences at high momentum transfer. For the magnetic form factor the 
results deviate only in the dip region. The dip position in the ET approximation is at too low momentum transfer 
as compared to the data. Since the dip region is sensitive to the MEC coupling constants, no specific conclusions 
can be drawn from this. In particular, adapting the unknown luej coupling constant one can readily reproduce the 
experimental data. Also the tensor polarization £20 are similar in both calculations. They agree with the experimental 
data (Fig. f7|). At q > 4fm~ 1 the predictions of the ET approximation are somewhat above the results of the BSLT 
approximation, which is due to a lower value of the dip in the charge form factor. 



VI. ELECTROMAGNETIC BREAKUP OF THE DEUTERON 



We now turn to describe the procedure to calculate the structure functions for inelastic electron scattering. There 
are two terms which are calculated separately. One is the PWIA and the other one is the FSI contribution. To 
determine the current matrix element we need to know the isospin structure. For the free final state we have both 
isospin I = or I = 1, whereas the deuteron is an isospin 7 = state. Clearly, the isospin has to be chosen for a 
given total spin S' of the two-nucleon system in accordance with the final state being totally antisymmetric. Since 
the outgoing nucleons are on-mass shell we have the conditions Pq — Ei P , +pf + Ei P ,_ pf and pf ■ P' = 0. In the 
following we consider the e.m. current matrix element with the final np pair in a total spin («S", M' s ) state. 



A. The PWIA contribution 



The simplest contribution to the deuteron current is the coupling of the photon to one of the two nucleons in the 
deuteron without FSI between the two outgoing nucleons. In Fig. || all possible contributions are depicted. Both 
initial and final states are antisymmetric and therefore the contributions (a) and (b) are identical. Let us now adopt 
the convention, where particle 1 i n the final state is assumed to be the proton, carrying momentum \P' + pf. To 
determine the structure functions ( 2.13 ) we calculate the e.m. current matrix elements in the cm. frame of the final 
state. In this frame we have Pq = 2E p f and pfo = 0. The PWIA contribution can be written as 



<P',p f ,S',M' s \J^ WIA \P,M> 



iey/2 

i—l 



x < P' ) p / ,5',M^|rW5W(pW,P)*g f) (pW ) P), 



(6.1) 



where = Pf — \q, and p^ — pj + \q. The state \P',pf,S',M' s > describes the free np pair with total spin S' 
and its explicit form is given in Appendix A. 

Since the deuteron vertex function is usually calculated in its rest system, it has to be boosted to the cm. system 
of the final np pair. The Lorentz transformation which does this is given in Eq. (2.S). The corresponding one for 
spin-i particles is given by 



A(£) = 



E D + M L 
2M d 



i + tV 



-Q 



E D + M D 



(6.2) 



Boosting the initial state to the cm. frame of the np pair, the deuteron current contribution Eq. (3.1) can be written 
in the form 



<P>, Pf ,S',M' s \J™ 1A \P,M> 



iV2 



E 



4ttM d 

i—l 

x < P',p f ,S',M' s \f®S®(k®,P)* L ?>(k< i >,P) i 



(6.3) 
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where f f = vfk{C), A(C) = A« (£)A( 2 ) (£) and JfcW = 

Introducing the deuteron state components ^ nm (p, p ) labeled by n = {J, M, L, S} and a combined helicity and 
p-spin label m = {A l7 A2, p}, the deuteron current Eq. ( |6.3| ) can now be evaluated in the same way as has been done 
in the elastic case. For a more detailed discussion we refer to jp|. On the chosen basis the deuteron current is found 
to be of the form 

<P\Pf,S\M' s \Jr iA \P,M>=^-± 

l—l 

x E E Dfr sy ,(n f )47^^ (6.4) 

m 1 pm.n 

where ml = {A'^A^H — h} and the explicit form of the one-particle propagator S^ p is given in Appendix A and 
Eq. (|]^). The momentum fcW is restricted by momentum conservation, since the spectator particle is on-shell we get 
that fcg 1 ^ = E ~ -^feC 1 ) anc ^ = — ^ + -^fe( 2 ) • I n wr iti n g down Eq. ( |6.4| ) we have shifted a 70 matrix to the vertex, 
so that the vertex operator is now given by f ft = rftA(C)"f^ and fft — T^A(£)7q 1 ^ . The extra 70 factor is due 
to the use of the closure relation for the helicity spinors. In Eq. (^J) the (^-dependence can explicitly be evaluated. 
For the spherical components (with p, = 0,±1) of the e.m. current we find that the total ^-dependence is given 
by exp[i(M + p — M' s )ip]. With this ^-dependence of the current operator, Rl and Rt don't depend on <p and is 
the ^-dependence of Rtl and Rtt given respectively by cost/? and 008 2^. The isospin structure contained in 

can readily be evaluated. Since particle 1 is assumed to be the proton we get rft -» F p /V2 and vft F n /V2. In 
the actual calculations it is convenient to treat the proton and neutron contributions to the e.m. current separately. 

Denoting the term where the photon interacts with the proton (neutron) as jj^ (j^^ we find with the help of 

symmetry relations for the deuteron vertex functions (see Q) 

J(")( P/ )H-1) L+S+5 '^ P) (-P/). (6.5) 

with all e.m. proton form factors in replaced by the corresponding ones of the neutron. 

The ET approximation for the PWIA contribution can be obtained immediately from the above expressions, by 
simply replacing the initial deuteron bound state by the BSLT vertex function. The relative energy variable in the 
propagator is prescribed by the condition that the final state describes on-shell particles. 



B. Rescattering contribution 



When FSI is included in the IA contributions, the structure is more co mp licated. Starting from the full Bethe- 
Salpeter formalism the current matrix elements has the form given in Eq. (4.4). We have 



<P', Pf ,S',M' s \J^ SI \P,M>= 



(6.6) 



Assuming again that particle 1 in the final state is the proton, we only h ave to take the two graphs (a) and (b) into 
account as represented in Fig. ||. Using the same analysis leading to Eq. (6^4) we find 



< P',p f ,S',M' s \J^ bI \P,M > 



— -^yy y y my 

x 1 i—l n' ,n m ,m' ,m" pp' 



(6.7) 



where the partial wave components 4'n';m',m"(PfiP') correspond to the isospin dependent t-matrix elements with total 
angular momentum n' — {</', M'} between the energy-helicity spin states m! = {X[, A' 2 , ++} and m" = {A", A2, p"} ■ 
The isospin dependence of T in Eq. (6/7) is for an isoscalar transition given by F s /y/2 and for an isovector transition 

(-iyF v /V2. 
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Replacing the integration variable p' by —p' in the integral we readily see that both contributions are equal. More- 
over, collecting all the (^—dependence we find e %lyM ~ M s)vf ^{m+^-m )<p ^ SQ fae (^'-integration can be performed 
with the result 

2ttS(M + h- M')e i( - M+ ^- M ^f' . (6.8) 

Notice that the same <p /-dependence exists as for the PWIA. In both the initial and final state the relative energy vari- 
able is set to zero, corresponding the BSLT prescription. Since the only j»Q-dependence is in the nucleon propagators, 
the ^-integration can be done analytical. As a result we are left with a two-dimensional integral which can be done 
numerically using standard gaussian quadratures. In so doing, one particular point has to be taken care of. At the 
point E p i = ^M np the particles are on mass-shell and we have a singularity in the positive energy propagator, which 
is of the form 1/(E P ' — \M np ). This singularity can simply be dealt with using a standard subtraction technique. 
Having dete rmined the e.m. current matrix elements in this way the deuteron st ruct ure functions can be constructed 



using Eq. (2J3). It should be noted that not all partial wave matrix elements in Eq. (6/7) need explicitly to be calculated. 



Due to parity conservation we have 

< J'M'L'S'\J^\P, M >= (-ly+J'+M+M'+n <J >_ m'L'S'IJ^P, ~ M > • ( 6 - 9 ) 

Besides this symmetry relation we have the selection rule M' = M + fi for the above matrix element, being a 
consequence of rotational symmetry. 



VII. RELATIVISTIC ANALYSIS OF INELASTIC ELECTRON SCATTERING 



In the preceding sections we have developed a relativistic framework to describe electrodisintcgration of the deuteron. 
The full relativistic form of the e.m. operator is employed and at the same time the relativistic structure of the nuclear 
dynamics is incorporated. To see what kind of predictions our fully relativistic formulation leads to in the case of 
inelastic electron scattering we study the structure functions for two kinematic situations. One is the inclusive 
experiment carried out by the Bates-Mit group and the other is the exclusive analysis by the Nikhcf-K group. Since 
we also want to compare our predictions to those obtained from nonrelativistic calculations, let us first discuss the 
e.m. operator used in the Schrodinger approach. 



A. Nonrelativistic limit 



In most of the studies of electron-deuteron scattering it is implicitly assumed that to a good approximation the 
deuteron can be described by a boundstate wave function satisfying a Schrodinger equation. In addition some effective 
nonrelativistic (NR) form is used for the e.m. operator. Various versions have been proposed in the literature based 
on taking the NR limit of the current operator. The e.m. operator is expanded in q/M, keeping only the first order 
term in the expansion and neglecting recoil effects. Because NR wave functions are used, boost effects are in general 
also neglected. The resulting effective e.m. operators have all in common to be the same in leading order, but to differ 
in relativistic order (q/M) 2 Since we have included in our study the full relativistic structure of the e.m. operator, it 
is obviously of interest to see how these relativistic corrections look like. In the calculations the vertex is defined in 
the cm-frame of the final state. Because the helicities in this frame and the cm-frame of the initial deuteron state are 
not equal, we have to replace the boost operators by Wigner rotations. A detailed discussion is given in Ref. ||. 

In the NR limit all particles are taken to be on-shell and the negative energy state components are neglected, i.e. 
since the nucleons are assumed to behave nonrelativistically, we may assume that the nucleons can be described by 
the positive energy spinors 



u{ 1 \p) = N p 



,(2), 



P) = N p 



1 

E p + M 
1 

Bp + M 



XX, 



X-A, 



(7.1) 



where N p — y/E p + Mjv /2E p and xx are the Pauli spinors with helicity A. Therefore the current matrix element to 
study in this approximation is given by 
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Ju 



ux>(p') (Fn M 



2M 



N 



-F 2 a^ v q v u\(p). 



(7.2) 



It can immediately be reduced to a Pauli form. We get 

1 



Jo = F 1 - ^t(Fi + 2F 2 )(q 2 + la ■ [K x q]) + O(^) 

(K + i[a x q])(l - -^{K 2 + q 2 }) - -L (q + i[a x K])(K • q) 



2M N 
F 2 



2M N 
1 

~8M 2 



i[cr x q](l 



8M : 



8M 2 



rq 2 ) 



2Miv 



8M 2 



(q + i[a x K]) 



([K x q] x q-i(cr-K)[K x q]) 



(7.3) 



with K = p + p' and where we have kept all second order terms. 

In the usual NR limit only the first order terms are kept in the charge and current operator. In so doing, we get 
the operator as used by Mathiot pa] and Sommer pq] 



Jo = F\, J 



1 



2M 



-FiK 



N 



2M 



N 



-G M (q x a). 



(7.4) 



Note that the correction terms in Eq. ([7j|) are a rather complicated function of both q and K. Except for the region 
of low q and small recoil corrections there is no a priori reason to expect them to be small. 

Another form which is often used, is found by starting from the alternative form of the relativistic current operator 
in terms of Sachs form factors 



G E = Fi 



4M£ 



F2, Gm — F\ + F 2 . 



Taking the NR limit in this case yields the e.m. operator used for example by Leidemann and Arenhovcl 



Jo — Ge 7 J — 



1 



2M 



-G F K 



N 



2M 



-G M (q x a). 



N 



(7.5) 



(7.6) 



The only difference between the two approximations is the use of F± or Ge in the charge and current op erator. From 



Eq. (7J3) we see that the difference is of order q 2 /M^, being of relativistic origin. Considering Eq. ( [7.3|) we see that 



the charge operator can be rewritten up to order q as 



Jo = G E (l 



SM 2 N 



1 



SM 2 N 



(2G M - G E )io ■ (K x q). 



(7.7) 



Momentum conservation at the vert ex g ives K = p + p' = 2p + q. Hence for small values of the missing momentum, 
i.e. p << q, the second term in Eq. (7/7) can be neglected. For values of q up to a few hundred MeV it suggests that 
the choice of Ge is indeed reasonable for the charge operator. The procedure as followed in Eq. (7.7) can not be used 
for the current operator, because of the lack of a small parameter. However, we can show that the isovector part up 
to order q A is given by 



2M N 



G M (1 



3q 2 

8M 2 , 



)(q x cr), 



(7.; 



which shows that the choice of Gm m Eq. (7.4) and Eq. (7.6) is reasonable for small q. It should be noted that in 



some calculations of Ref. 



the factor G, 



Up 



EnEj 



has been added to the cross section, whereas this term is already 



included in Eq. (7.2). This correction is clearly of relativistic order. 

Besides the additional terms neglected for a given choice of the NR e.m. operator, other relativistic effects like 
contributions from the boost transformation and the negative energy state components of the deuteron and scattering 
wave function, can play an important role. In Fig. ^ are shown the relative contributions from these corrections to the 
response functions for an exclusive experiment with missing momentum p m = 100 MeV/c and q in the range of 300 
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MeV/c to 1500 MeV/c. In addition we used 6 P — 120° and 9 e — 60° and took the kinematics in plane. In this figure 
we are comparing the different approximations to the complete relativistic calculation with FSI included, but without 
negative energy state contributions. The effect of including the negative energy state contributions is given by the 
double dot-dashed curve. As a second approximation we may drop the boost operators in the current op erat or (solid 
line). This is similar to the NR approximations except that we keep the complete current operator Eq. (7.2). If only 
the PWIA contribution is considered all response functions show identical dependence on the boost operators. The 
contributions from the boost to Rl,Rt and Rtl become more important at higher q values and can be as large as 
40 % around q — 1500 MeV/c, whereas in the case of Rtt inclusion of the FSI substantially enlarges this dependence 
and the boost effects are large even for low momentum transfer. The negative energy state contributions increase at 
higher q and are of the order of 5% around q = 300 MeV/c. Again a different behavior is found for Rtt, an effect 
of around 30% is seen for the whole kinematic region. An important conclusion is that the negative contributions 
cancel to a large extent the boost effects. Similar as for elastic electron deuteron scattering a consistent relativistic 
treatment is also important in this case. 

Next we may consider the effect of neglecting the FSI contributions. This yields the long dashed curve in Fig. || As 
expected we see an increase in the FSI contribution at lower q, since q is directly related to the energy of the np-pair in 
the final state. Also it is of some interest to see the sensitivity on the choice of the relativistic two-nucleon propagator 
on the response functions. The dashed line depicts the relative difference between the BSLT and ET approximation, 
including FSI. Except for Rtt, the difference is less than 5% and increases for higher values of q, whereas for Rtt it 
decreases with q. This dependence on the propagator structure only exists if FSI is inclu ded. 

Finally, also plotted in Fig. ||are the predictions with the e.m. operators, given by Eq. ( |7.4D and Eq. ( |7.6| ) with the 
Fi and Ge form factors. The results are respectively given by the dotted and the dot-dashed lines. As expected we 
see increasing deviations with increasing q. For Rt the NR Ge operator yields clearly the closest resem blan ce to the 
calculation using the complete current operator. The latter is in accordance with our discussion of Eq. ( |7.7| ) that the 
use of Ge in the charge operator should be reliable in the quasielastic region. Beause of isovector dominance in the 
contribution of Rt both NR operators give the same results and deviate only a few percent of the results with the 
exact current operator, what is in accordance with Eq. ( [7. 8[ ). Finally we see from the figure that the deviations are 
very large for the response functions Rtl and Rtt even at comparatively low momentum transfer. These deviations 
are essentially due the interference structure of the response functions, therefore probing those parts of the current 
operators not described well by the NR operators. 

where we have kept all second 



Not shown in these figures are the results with the current operator of Eq. (7.3 



order terms. Similar higher order corrections have also been considered recently [28 29 For Rt and Rt although 
the results agree reasonably well with the full calculations at moderate momentum transfer they are found to be 
substantially different at higher momenta. In general at high momentum transfers the deviations are such that the 
complete e.m. operator including boost effects is needed for a reliable analysis of the response functions. For Rtt 
and Rtl calculations with this second order operator show it gives comparable results as with the complete operator. 
Deviations are within a few percent. At these low missing momenta no large model dependence of the wave functions 
is expected. The deuteron wave function is mainly S-wave in this region. We have tested the aspect of the non 
relativistic reduction of the propagators. With this approximation we find effects smaller than 0.5 %. 

get some insight about sensitivities on FSI and negative energy states in a larger part of the phase space in 
J and nOlcontour plots are shown where we have varied in addition 9 p . For the relative contribution of the FSI, 
shown in FigT^j, we find very similar effects in the various kinematic regions for Rt and Rt- From these figures we 
see that both contributions to in particular the Rtt response function can be substantial. 



B. Structure functions of the exclusive reaction 



As an example we discuss the exclusive experiments carried out at Nikhef-K |lj|2|] . In these experiments in addition 
to the measured final electron also the outcoming proton is d etec ted. As a result the cross section is described by the 



four response functions Rl, Rt, Rtl and Rtt given in Eq. (2.6). Rtl and Rtt have an out of plane dependence of 
respectively coscj) and cos2(f>. Consequently Rtl can be separated by carrying out two measurements at (f> — and 
4> = 7r, such that p • q is constant. The response function Rtl is interesting because through its interference structure 
it is sensitive to the isovector part of the charge operator and the isoscalar part of the current operator. By measuring 
only Rl and Rt no information is obtained about these terms because of the dominance of the isoscalar charge and 
isovector current operators. The experiments were performed with a missing momentum p m in the range of 40 to 
180 MeV/c. In this kinematic region the FSI is expected to give only small contributions due to the low missing 
momentum and therefore the high relative momenta of the outcoming nucleon pair. Because of the dominance of 
the PWIA contribution the chosen kinematics is very suitable to examine the underlying structure of the one-body 
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current operator. 

In Fig. [n] the results of the Nikhef experiments together with our predictions are shown. In all these calculations 
the OBE model has been used to construct the continuum and deuteron vertex functions. In these figures PWIA 
corresponds to our relativistic PWIA calculation and the curve denoted by REL(-) contains in addition the FSI (both 
positive and negative energy states). For both Rl and Rt the FSI is small and increases with increasing missing 
momentum. Also are plotted the predictions obtained using a nonrelativistic description with e.m. operators, given 



by Eq. (7.4) and Eq. ( |7.6| ) with F\ and Ge form factors. The results show clearly that the current operator with 
F\ is too high as compared to the experimental data for Rl and the Ge description gives a much better agreement 
with both the data and the relativistic calculation. The results for Rtl and Rtt are remarkable. Although we are 
at relatively low momentum transfer, relativistic effects are substantial. We find that our relativistic calculations 
leads to a good description of the Rtl data, whereas both NR predictions are clearly too low. In this case the NR 
F\ operator is closer to the relativistic prediction. Similar conclusions can be drawn for the out of plane structure 
function Rtt- It is worth noting that for this case the FSI effect is considerably larger. 

From these results we may conclude that in the kinematic region considered the most important relativistic effects 
arise from the higher order corrections to the e.m. operators employed in the usual NR analysis. In general the NR 
descriptions are clearly not adequate and they can represent a poor approximation to the full relativistic operator. 

C. Inclusive scattering in the quasielastic region 

We now turn to describe our calculations for the inclusive reactions in the quasielastic region. These experiments 
were done at Mit-Bates for momentum transfers q in the range of 300 — 500 MeV/c 0. The quasielastic peak is 
defined by 



M D +w= ^M 2 n + q 2 + M N . (7.9) 

In the considered exp eriments a longitudinal-transverse separation has been carried out. The four momentum transfer 
is such that Eq. (73) to a good approximation reduces to to = q 2 jlM^. For these values of u> and q the electron is 



essentially probing the nucleons with missing momentum near zero. 

In Figs. |l2f]l4| are shown the various corrections to the longitudinal and transverse response functions at momentum 
transfer q = 300, 400 and 500 MeV/c. The full relativistic expression of the electromagnetic vertex is used. Inclusion 
of the boost operator shows only a very small contribution. Considering Rl we see that the PWIA describes to 
a good approximation the response, except at lower momentum transfer. The Born contribution can essentially be 
neglected, because F" is very small. On the other hand the FSI is seen to be important for Rt over the whole range of 
considered momentum transfer. Moreover, since the magnetic form factor of the nucleon is comparable in magnitude 
to F% the Born term also gives a substantial contribution to Rt- For larger ui the delta-degree of freedom is seen to 
come up. This contribution and the other MEC's will be discussed in the same relativistic framework in a next paper. 

Including only respectively the Born, the J = and 1 NN scattering amplitude in the two response functions we see 
that most of the FSI has been accounted for . The main FSI contribution to Rl is from the J — 1 ( 1 Pi, 3 Si~ 3 D\, 3 Pi) 
amplitudes. Especially at low ui this contribution can be significant due to the decreasing relative momentum of the 
outcoming nucleon pair. Similarly as for the three-nucleon systems [30| l the response increases at low u>, whereas it 
decreases at the quasielastic peak due to the FSI contribution. The PWIA result for Rt is clearly below the data. In 
particular, the J — ( 1 S' , 3 P ) FSI contribution is only of some importance for Rt- Due to the FSI contributions 
our overall predictions are well in accordance with the experimental data. Although substantially smaller than in the 
the trinucleon system some strength tends to be missing in the transverse response function at lower energy transfer. 
We also considered the negative energy state contributions. Exactly as we have found for exclusive scattering this 
contribution seems to be vanishing at the quasielastic peak. Furthermore it reduces the strength at low u>, whereas it 
tends to increase slightly the response functions past the quasielastic peak. 

The overall agreement with the experimental data is gratifying. Since at these relatively low momentum transfers 
the existing parameterizations of the e.m. nucleon form factors are very similar, less dependence is found on the 
choice of these nucleon form factors. Using for example the Iachello et al |3l[] form factors changes of only at most 1 
percent were found. 

VIII. SUMMARY AND OUTLOOK 

We have presented in this paper a framework to study the e.m. deuteron current based on a relativistic quasipotential 
approach. The nuclear dynamics itself is described in terms of a symmetric choice for the relative energy variable. 
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Assuming that the two-nucleon vertex functions depend smoothly on this variable, the current matrix element is 
determined by explicitly integrating out the relative energy variable in the e.m. vertex loop. This equal time 
approximation has the advantage that it can be used to describe both elastic and inelastic electron scattering and 
that it in principle can be extended in a systematic way. As a direct application the elastic e.m. properties of 
the deuteron are studied using this formalism. The differences between this approach and the one used in Ref. Q 
for elastic electron deuteron scattering are the treatment of the initial and final states and the nucleon-propagator 
structure of the one-body current. We have shown that for elastic scattering these two approaches give compatible 
results and differences are found to be minor. 

The approach discussed here has clearly the advantage that besides elastic electron scattering case, it can also be 
applied to the deuteron breakup in electro- and photo-disintegration processes. Explicit expressions are written down 
for the deuteron current in the case of inelastic electron scattering. The helicity formalism is extensively used to relate 
the scattering wave functions of the relativistic one boson exchange model to the current matrix element. Since there 
are also isovector current contributions in this case, mesonic exchange currents are in general needed to satisfy gauge 
invariance. These currents can also be constructed in this ET approximation and will be studied in a subsequent 
paper. 

Various different forms of one-body operators exist in the literature to describe the e.m. processes in a nonrelativistic 
framework. These are derived from relativistic operators which are on-shell equivalent, but lead to different operators 
in the NR limit. We have compared the NR calculations with our relativistic description of the electro disintegration 
of the deuteron. The important difference is that we also take in our approach the higher order terms in q/M into 
account. It should be noted that we have assumed a specific on shell form of the e.m. nucleon form factor and not 
considered its off-shell dependence. For this a detailed dynamical model of the nucleon is needed. 

One interesting experiment to study the difference between a relativistic and nonrelativistic analysis is the exclusive 
experiment performed at Nikhef-K, in which the three response functions Rl, Rt and Rtl could be separated. Due to 
the small missing momentum in these experiments the FSI-contribution is very small and the response functions are in 
a good approximation described by the one-body current (PWIA) . In our analysis we have found that both relativistic 
and nonrelativistic theories give a good description of Rt- Differences show up for Rl, the NR-theory using the Ge 
form factor gives a good description, while the F± current operator is systematic to high. Interesting results were 
found for the interference structure function Rtl- Both NR-approaches failed to describe the data. Including the 
next order in q/M to the one-body current operator showed a much better agreement. Other relativistic effects have 
also been considered. These are due to boost effects and the presence of negative energy states. In the kinematics 
studied these contributions were found to be small. However, it is expected that these corrections will be substantial 
at higher momentum and energy transfer as can be seen in the case of elastic electron scattering. We have tested 
this assumption by considering the range of q = 300 MeV/c to 1500 MeV/c. Effects of 40% were found. The second 
experiment we discussed in our relativistic approach is the inclusive experiment done at Bates. In this experiment 
much larger effects due to the Born term and FSI are found. A good agreement with the data is found, although at 
larger momentum transfer q the longitudinal response function is slightly overestimated. This may however readily be 
explained in terms of model dependence on the deuteron structure. For these experiments the NR-description using 
Ge gives similar results as the relativistic calculation. 



APPENDIX A: HELICITY BASIS 



In this Appendix we summarize the helicity basis used in the partial wave decomposition of the NN amplitude. For 
more details we refer to Refs. Following Kubis @ we use for the positive and negative energy spinors 



u[ + \p) = N p 
u{-\p)=N p 



1 

2\p 

E p +M 

-2\p 
E p +M 
1 



(Al) 



where xa($, f) are the two-component Pauli spinors, which are given by 



cose/2 \ 
e lv sin 6>/2 ) ' 

-e - ^ sin 61/2 
cos6>/2 



(A2) 
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These spinors are normalized as 

4 p)t (p)4^(p)=V^A<, (A3) 



where p — ± and the normalization is given by N p — \/{E p + Mn)/2E p . Let us consider the projection operators 

A P (p) = E u i P) (P)4 rt (p)- (A4) 

A 

We may rewrite the single nucleon propagator 

S(p) = ^-M N + ier = /j^f + . e (A5) 

in terms of these projection operators. Using the relation 

i> + M N = (E p +j3 )A+(p) + (po - £ ; P )A-(p) (A6) 

we get for the propagator 

S(p)= A+(P) + A - (P) . (A7) 
Po - E p + ie po + E p + ie 

We now turn to discuss briefly the partial wave representation of the two- nucleon t-matrix. When two-particle states 
are involved, we take the convention of Jacob and Wick Q where particle 1 is described by the spinors u x P \p), while 

for particle 2 the spinors U_l(— p) are used. The two-particle states built from these helicity spinors form a complete 
basis \pX1X2p > (with p = (pi, pz)) in Dirac space. In this representation the t-matrix is given by 

<p f X[X' 2 p'\ct>(pf,p,P)\pX 1 X 2 p> . (A8) 

Its angular dependence can be exhibited using the total angular momentum states 



I JMA1A2 >= / dnpD{; x (9,ip)\p,6,i P ,X 1 ,X 2 > . (A9) 



In view of rotational invariance of the t-matrix, Eq. (A.8) can be written in the form 

2 7+1 

t = J2<t>n;m',m(pf,p) = £ — ^^(^^.mCP/.P^MA^). ( A 10) 
n J,M 

where we have defined 4>m',m(Pf>P) =< PfJ^i^P'l't'iPf^P^^pJ^i^P > with n — {J,M},m' = {A'^A^p'} and 
m = {Ai, X-2,p}- These angular momentum states can be related to states labeled by the set {J MLS}. We have 



J,M L'.S'.L,S 

xDi' IX ,(n f )C^; J Di; x (Q p )C^ s x J < JL'S'\tf> J (j>f,p, P)\JLS > . (All) 



A detailed analysis is given in the Refs. |6|, |14| , |32| . The free np pair can be described by the two-particle states |pAiA2 > 
with p = (+, +). These helicity states are connected to the total spin S states |p, S, Ms > through the relation 

|p, S, M S >= ^£a(^ p )C a Y_ 5 A2A |pA 1 A 2 > . (A12) 

Ai A 2 
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APPENDIX B: THE K LOOP INTEGRATION IN THE E.M. CURRENT 

In the current operator for the IA in the equal time approximation, initial and final state are assumed to be 
independent of the relative energy variable. As a result the only dependence on this variable is in the nucleon 
propagators. Let us first consider the inelastic case. The integral we want to study is of the form 

I = Jdk>S 2 (k',P>)f»S^(k,P) (Bl) 

where T M is the e.m. operator including the boost operators. Both sets of variables (k 1 , P') and (k, P) are in their 
own cm. frame, i.e. P' = (0, 2E') and P = (0, 2E) where 2E' = M np and 2E = Mp. The Lorentz transformation C 
of the cm. frame to the lab frame is given by Eq. ( |2.9| ). 

To perform the k' integration we have to analyze the position of the singularities in the propagators. The singu- 
larities in the two-particle propagator S2(k,P') are located at 

(l)k' = E' + E k , - is, (2)k' = -E' + E k , - ie, 
(3)k' =E'- Ev + ie, (A)k' = -E' - E* + ie, 

(B2) 

while in the initial one-particle propagator S^(k,P cm ), expressed in the final state variable k' , are given by 

(6)k' = {oj cm - E 1 ) - ££"_ q + ie. (B3) 

1(2) /(3) 

In the breakup region E > a pinching can occur between the poles fc and fc . This generates the elastic 

cut defined by E^ = E' . Moreover, the possibility exists that the singularities k'^ and k'^ may cross. This occurs 
when 

- E k , = -uj cm . (B4) 

In Fig. [l5| is shown the kinematic region of q and uj where this can happen. The solid line is for small u> given 
by uj w q lab /2Mm which corresponds to quasielastic scattering. Note that for real photon absorption uj = q this 



double pole doesn't occur. The condition (B4) corresponds to the situation when the hit particle before and after the 
absorption of the photon is on the mass shell. Since these poles are on the same side of the contour, their crossing 
doesnot correspond to a pinching singularity. 

By closing the contour in for instance the lower half plane the integral over k' can be performed. The integral picks 

up contributions from the poles at k' , k'^ and k'^ . The singularity at kg ^ is a negative energy state contribution, 
the other two are from positive energy states. In so doing we get 

/ = 27™A 2 -(k')S (2) (fc'W, P')f At 5 (1) (fc (1) , P) 
-27nA+(k')S^(k'^,P')r^S (1 \k^,P) 

-2mSM (k' ( V , P')S& (fc'(5) , P ')f , A+ (k) -H- , (B5) 



TPC 



where the four vectors k'^> and k^ have as fourth component k^ and kg respectively. If we consider only positive 



energy intermediate states Eq. (Bl) reduces to 

7=-2 7 riA+(k')A+(k')f /J A+(k) 

1 ' 1 1 ^ , (B6) 



E + fc<, 2) - E k m - E k , ) E , + _ Ev E , _ k f) _ Ek , E k , _ q J 



It should be noted, that both contributions develop separately a double pole for k satisfying Eq. (B4). However, the 
two terms cancel at this point leading to only a single order pole. Furthermore notice that only the first term has a 
pole of the form (E' — Ek')" 1 which generates the elastic cut. 

In elastic scattering we take as integrat ion var iable ka and the the vertex is evaluated in the Breit frame. In 



terms of the variable ko the singularities (B2IB2I) are given by 
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(1) fco = E + E u -ie, 

(2) fc = -uj - E + £ k + q - ie, 

(3) fc = £ - # k + ie, 

(4) fe = -u-E- -E k +q + «£, 

(5) fc = -E + E u -is, 

(6) fc = -£ - £ k + «£• (B7) 

It can readily be shown that the kinematics in this case is such that the poles k^ and cannot cross. Closing the 
contour in the lower half plane, the integral becomes 

I elastic = 2wiS^(k'^,P')t^ 1 \k^,P)A^(k) 
-2mS {1) (fc ,(5) , P')t „A+ (k)5 (2) (fc (5) , P) 

-2ni-^A+(k')f^ 1 \k < - 2 \P)S^(k ( - 2 \P). (B8) 
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FIG. 1. The phase shifts for the BSLT equation as a function of the lab energy for the fits given in Table I. The solid lines 
are the calculated phase shifts for fit B with the negative energy intermediate states included. The dashed and dotted-dashed 
lines correspond to fits A from Ref. |J and B respectively where only positive energy intermediate states have been kept. The 
data are from Ref. fllq |. 

FIG. 2. The PWIA (a and b) and the Born (c and d) contributions to the deuteron current. 

FIG. 3. Feynman diagrams corresponding to the IA contribution to the deuteron current with FSI included. 

FIG. 4. The charge form factor Fc and the quadrupole form factor Fq in the relativistic impulse approximation using the 
ET and the BSLT prescription for the deuteron current. The data are from Ref. po| . 

FIG. 5. The same as Fig. ^, but for the magnetic form factor Fm- 

FIG. 6. The full predictions using the ET and the BSLT prescriptions for the deuteron current with and without MEC 
contributions, but for the electric and magnetic form factors A and B. The data are from Refs. J2lj and p2| . 

FIG. 7. The same as Fig. ^, but for the tensor polarization t2o- The data are from Refs. [ po| , p3|j24| ] . 

R a (r ) r (r ) 

FIG. 8. The relative difference q oixr (r ") °^ ^ ne res P onse functions R a without negative energy state contributions 
and R2, calculated in various approximations. When R a is the full result we find the solid curve. Evaluating R a by neglecting 
both the negative energy states and the boost transformations yields the double dot-dash line. Dropping furthermore also the 
FSI contributions in R a yields the long dashed curve. The effect of the BSLT propagator choice is given by the dashe d cu rve. 
The dotted and dot-dashed curves correspond to the results of the relative change for the various NR approximations 7.4 and 



7.6 to the current. 



FIG. 9. The relative contribution of FSI to the four response functions as function of q and cos(6 l ). 
FIG. 10. The relative contribution of negative energy states to the four response functions as function of q and cos(#). 
FIG. 11. The response functions for the exclusive Nikhef kinematics. The data are from Ref. M. 
FIG. 12. Longitudinal and transverse response functions for q = 300 MeV/c. The data are from Ref. ||. 
FIG. 13. Longitudinal and transverse response functions for q = 400 MeV/c. The data are from Ref. ||. 
FIG. 14. Longitudinal and transverse response functions for q = 500 MeV/c. The data are from Ref. JBJ. 

FIG. 15. In the region I there is a crossing between two positive energy poles of the one-particle propagator were before and 
after absorption of the photon the particle is on mass shell. In region II this crossing doesn't occur. 
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